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ABSTRACT 

All far ultraviolet observations of HD209458 tend to support a scenario in which the inflated hy- 
drogen atmosphere of its planetary companion strongly absorbs the stellar Ly-a flux during transit. 
However, it was not clear how the transit absorption depends on the selected wavelength range in the 
stellar line profile, nor how the atomic hydrogen cloud was distributed spatially around HD209458b. 
Here we report a sensitivity study of observed time and spectral variations of the stellar flux. In 
particular, the sensitivity of the absorption depth during transit to the assumed spectral range in the 
stellar line proflle is shown to be very weak, leading to a transit depth in the range (8.4 — 8. 9)% ±2.0% 
for all possible wavelength ranges, and thereby confirming our initially-reported absorption rate. Tak- 
ing the ratio of the line profile during transit to the unperturbed line profile, we also show that the 
spectral signature of the absorption by the exoplanetary hydrogen nebula is symmetric and typical of 
a Lorentzian, optically thick medium. Our results question the adequacy of models that require a huge 
absorption and/or a strong asymmetry between the blue and red side of the absorption line during 
transit as no such features could be detected in the HST FUV absorption profile. Finally, we show 
that standard atmospheric models of HD209458b provide a good fit to the observed absorption profile 
during transit. Other hybrid models assuming a standard model with a thin layer of superthermal 
hydrogen on top remain possible. 

Subject headings: Stars: individual (HD209458) — Stars: planetary systems — Ultraviolet: stars — 
Line: profiles — Techniques: spectroscopic 
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1. INTRODUCTION 

Just like our Sun, the Ly-a line profile of HD209458 
shows a double-horn emission with the distinguishing 
feature of intervening interstellar gas strongly absorbing 
the line c ore. From t he signal left at the Earth's orbit, 
[Vidal-Mad Tar et al.. I (]2003) (hereafter VM03) initially 
claimed that HST/STIS medium resolution observations 
of HD209458 are showing a strong absorption of 15% 
during transit and that a large asymmetry appears on 
the blue side of the Ly-a absorption profile. Those 
authors used the 15% drop-off of the stellar signal during 
transit as an indication of the presence of a hydrogen 
cloud that extends beyond the Roche lobe of the planet. 
The next step for VM03 was to use the asymmetry 
claimed in the blue side of the line during transit as an 
indication that the evaporated gas is accelerated by the 
radiation pressure from the star, an acceleration that 
results into the build-up of a cometary shape of the 
nebulosity that escapes from the planet away from the 
star. 



Later on, iHolmstrom et aO ()2008f l rejected the 
acceleration process by radiation pressure, arguing that 
neutrals will not have the time to be accelerated by the 
radiation pressure before they get ionized. Furthermore, 
they proposed that observations thus far reported by 
VM03 on HD209458b could be explained by the presence 
of a population of energetic neutral atoms produced 
by charge exchange between the neutral atmosphere 

^ Visiting Scholar, Department of Applied Science, UC Davis, 
CA 95616 



and the stellar wind ions beyond an obstacle region 
corresponding to the planetary magnetopause assumed 
around ~ 4.2 Rp from the center of the planet. In the 
meantime, we revisited the HST/STIS medium resolu- 
tion archive data to report that exoplanet HD209458b 
absorbs ^ 8.9% ± 2.1% of the parent star Ly-a fiux 
during its tra nsit, mu c h less than previously reported 
by VM03 ( Ben-Jaffel. I (l200l . hereafter BJ07). We also 
reported that no sign of asymmetry could be found in the 
stellar line profile during transit, a finding that directly 
contradicts the claims initially reported in VM03. Soon 
thereafter, the absorption drop-off reported by BJ07 for 
a specific wavelength range was used out of the spectral 
information provided in Fig. 3 of BJ07 to claim that 
evaporation was strengthened with a factor ~ 2 asym- 
metry in the absorption between the blue (9.8% ± 1.8%) 
and red (~ 5.2% ± 1.0%) sid es of the Ly-a line during 
transit (jVidal-Madiar et al. I (|2008) hereafter VM08). 
In any case, it appears that our initial selection of the 
wavelength window to estimate the transit absorption 
created some confusion as to the reading of the data 
analysis presented in BJ07, and this confusion adversely 
affects efforts to discriminate between the processes thus 
far proposed to explai n the hyd r ogen d i stribution around 
HD20 9 458b ('VM03. IMuiioz. I (I200I : ISchneiter et aO 
(|2007f) : IHolmstrom et al.. I (|2008f )). In the following, we 
revisit available Ly-a HST/STIS medium resolution 
archive data to conduct a sensitivity study that should 
dissipate the existing confusion, an important step 
in formulating clear observational constraints on the 
steady-state properties of the hydrogen nebula around 
HD209458b. 
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Data analysis is carefully described in section 2 so that 
light curves may be extracted (in selected wavelength 
windows within the stellar line) and studied in section 
3. In section 4, time variability of the HD209458 Ly-a 
emission is examined and further evidence is provided on 
the variable behaviour of the HD209458 emission line, 
a property that may not be well-described by proxies 
derived for the solar Ly-a line. The absorption spectral 
profile is then derived as the ratio between the stellar flux 
during transit to the unperturbed spectrum of the star 
with no transit effect. Finally, the adequacy of models 
thus far proposed to explain the hydrogen distribution 
around hot giant exoplanets is discussed in the light of 
the constraints derived here, with new solutions proposed 
for the H nebula distribution. 

2. OBSERVATIONS AND DATA ANALYSIS 

Four HST visits during the planetary transit of 
HD209458 arc considered, each composed of two to 
three long exposures (^ 2000 s) distributed in time over 
the transit period (e.g.. Table 1). All exposures were 
obtained with the HST/STIS G140M medium resolution 
grating using the 52"x0.1" slit in the time-tag mode. 
It is important to emphasize that the time-tag mode 
is a powerful tool that keeps track every 125 x 10~^ s 
of photon events during each exposure. This mode is 
generally used to improve the spectral and imaging 
resolutions by correcting for the different moti o ns th at 
may affect the recorded signal (jBrown et al.. I (|2002f )). 
Doppler shifts due to HST motion, thermal drifts, and 
spacecraft jitter/drift are the main sources of motions. 
In the standard mode (ACCUM) used in several studies 
(VM03, VM08), correction is made only for the Doppler 
motion. For the other sources of motion, we used the 
time-tag mode, which offers the possibility of cancelling 
these motions by taking short sub-exposures that may 
be cross-correlated. In addition, with its high time 
resolution, the time-tag mode also offers the opportunity 
to monitor the time variation of the recorded signal, 
allowing the study of fast events and/or accurate 
correction for variable features by also considering short 
time exoposures. 

Following BJ07, our set of eleven exposures are 
sampled here using a time bin of 300 s from the avail- 
able time-tag information. Next, each sub-exposure 
is calibrated through the STIS pipeline. Because the 
HST/STIS/G140M grating is a spectro-imaging instru- 
ment, every exposure provides an image of the sky along 
the slit and the corresponding spectra on the cross dis- 
persion direc tion (Fig. 1). Comp ared to high resolution 
observations (|Wood et al.. I (|2005f )). the overlap between 
the stellar and sky signals is worse in the medium 
resolution, thereby making the corresponding correction 
more difficult (Fig. 2a). On average, the altitude of 
the HST orbit is high (- 500 Km) but stiU within the 
Earth's hydrogen geocorona. Therefore, any observation 
is contaminated by the emission and absorption from 
that upper layers of the Earth's atmosphere and by the 
absorption and emission from the interstellar medium. 
Insofar as the star is a point-like source that has few 
pixels extented over the slit spatial direction (Fig. 1), 
the other areas of the slit have only the sky background 
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Fig. 1. — A sample 2D Ly-a spectro-image of HD209458 as ob- 
tained by HST/STISG140M medium resolution grating using the 
long 52"x0.1" slit. The STIS MAMA detector has a ~ 25" X 25" 
size, making each pixel covering 0.024" over the sky. The dispersion 
direction is horizontal so each line of the image is a spectrum. The 
spatial direction is vertical. The position of the stellar spectrum 
along the slit is indicated by an arrow. The bright vertical band 
is the sky background distribution over the sky. Two ga ps, corre- 
sponding to the STIS fiducial bars IIBrown et al.. I l|2002l) ). appear 
at the bottom and top of the slit. 

emission which may a priori be used to correct for that 
contamination. As shown in Figure 2b, the difhculty 
of this correction comes from the uncertainty about 
whether the geocorona's signal strength when estimated 
from a detector sector along the STIS slit is different 
from the one where the stellar signal was recorded. 
Comparing the sky background signal from different 
sectors of the detector along the slit and for different 
conditions of observation, we confirm the conclusion of 
BJ07 that the STIS MAMMA detector has an inherent 
non-uniformity corresponding to an incompressible un- 
certainty of 5% on extended sources that is comparable 
to photon statistical errors. 

At this level, the correction for the sky contamination 
is not yet complete. Indeed, we still need to correct for 
the absorption by the geocoronal gas. While it may not 
be necessary when the target's signal is nearly zero in 
the wavelength range around the spectral position of 
the Ly-a line in the reference frame tied to Earth even 
so, we must check how the line spread function (LSF) 
of the instrument expands the absorption signature over 
the spectral domain of interest. In addition, the spectral 
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Dataset name 


Program ID 


Date Obs. 


Time Obs. 


Start time - TCT (s) 


Duration (s) 


Earth's Velocity 


O4ZEA4010 


7508 


2001-08-28 


06:50:07 




-7980.28 


2600. 


2.05 


O4ZEA4020 


7508 


2001-08-28 


08:26:25 




-2202.29 


2600 


2.05 


O6E201010 


9064 


2001-09-07 


20:35:27 




-8075.11 


1780. 


-2.59 


O6E201020 


9064 


2001-09-07 


22:12:36 




-2246.12 


2100. 


-2.62 


O6E201030 


9064 


2001-09-07 


23:48:54 




3531.89 


2100. 


-2.65 


O6E202010 


9064 


2001-09-14 


21:11:51 




-10167.67 


1780. 


-5.68 


O6E202020 


9064 


2001-09-14 


22:42:08 




-4750.68 


2100. 


-5.71 


O6E202030 


9064 


2001-09-15 


00:18:24 




1025.31 


2100. 


-5.74 


O6E203010* 


9064 


2001-10-20 


02:50:03 




-11258.45 


1780. 


-19.30 


O6E203020* 


9064 


2001-10-20 


04:19:45 




-5876.47 


2100. 


-19.32 


O6E203030* 


9064 


2001-10-20 


05:55:59 




-102.45 


2100. 


-19.34 



TABLE 1 

HST/STIS DATA SET ON HD209458 used in this study. All observations were obtained with the G140M grating and the 

52"x0.1" LONG slit. Transit central time (TCT) is defined by 2,452,826.628521 HJD iIKnutson et al.. I (120071) '). Earth's 

velocity is projected along line of sight in the heliocentric reference frame. (*) indicate data set that shows most time 

variability. 



position of the geocoronal absorption changes from 
one exposure to another in the hehocentric reference 
frame depending on the day of observation (see Table. 
1). Assuming standard models for the Earth's upper 
atmosphere, the Ly-a geocoronal opacity was evaluated 
for every exposure before the corresponding absorption 
profile is derived, convolved with the STIS/G140M LSF, 
and then applied to the sign al after the sky bac k groun d 
emission was subtracte d (iBen-Jaffel et al.. I (|2007t ): 
IChamberlain fc Hunten. I ()1987f) V Fortunately, most 
of the geocoronal absorption takes place within the 
spectral domain that we excluded from our data analysis 
so that its impact is less than 1% in the useful domain. 
Finally, because the interstellar absorption does not 
really change from one exposure to another, and should 
not affect our analysis of the transit effect, no correction 
is applied for this feature.. 

After all corrections were applied, the whole data set 
resulted in 53 bins time series of the HD208458 planet- 
star system versus the orbital phase angle that is quite 
similar to our initial time series within the range of sta- 
tistical errors (e.g., Fig.[l]of BJ07). To ensure that all the 
effects cumulated (LSF, long slit, thermal motion, etc.) 
did not corrupt the stellar signal per wavelength pixel at 
the 1% level, the wavelength domain of contamination 
by the sky background was defined as [121.541, 121.584] 
nm, the same as in BJ07. 

3. LIGHT CURVES SENSITIVITY TO WAVELENGTHS 
DOMAIN 

To obtain the trend of the planetary transit in Ly-a , 
BJ07 integrated the time series spectra in the range 
[121.483,121.536] nm on the blue side of the stellar 
line together with the range [121.589, 121.643] nm 
on the red side. Accumulating the signal inside the 
planetary transit period, BJ07 previously derived a 
^ 8.9% ± 2.1% drop-off of the stellar Ly-a intensity. 
The question then is: does the absorption rate depend 
on the selected wavelength range? For example, using 
low resolution HST / STIS observations of HD209458, 
IVidal-Madjar et al. I (|2004f ) found an absorption rate of 
5.2% ± 2 % for neutral hydrogen, a rate much smaller 
than they derived when using HST/STIS medium 
resolution observation. Recently, VM08 claimed that 
the difference between our estimate (8.9% ± 2.1%) and 
their initial absorption rate (15% ± 4%) comes from 
the distinct integrated wavelength ranges considered. 



These authors also claimed that they can reconcile 
medium- and low-resolution observations using the 
same argument of difference in the wavelength range 
namely, full line in the low-resolution case compared 
to restricted windows around the stellar line peaks in 
the medium-resolution observations analysis. Therefore, 
it is important to verify these claims and to estimate 
the sensitivity of the transit absorption to the selected 
wavelength range. 

Starting from the time series of the Ly-a line emission 
(e.g.. Figure 1 in BJ07), we have generated light curves 
(LCs) for several wavelength windows, including those 
considered by VM03 and more recently by VM08. As 
an initial step, we tried to compare LCs corresponding 
separately to the blue and red line sides (Fig. 3). 
Accumulating the signal from three bins inside transit 
with the corresponding time weights, our analysis shows 
an absorption ~ 8.5% ± 2.5% for the blue side and 
~ 8.2% ± 2.5% for the red side (see Table 2 for details). 
As shown in Figure 3, the two LCs behave similarly over 
time, except for the last bin around time ~ 5000 s where 
the red side seems even more absorbed. Despite the 
noise level, the observed trend for the two LCs supports 
the idea that both sides of the stellar line are absorbed 
almost equally during transit. By contrast, for the same 
wavelength domains, VM08 recently reported absorp- 
tions of 9.8% ± 1.8% and - 5.2% ± 1.0%, respectively, 
for the blue and red sides of the line, suggesting a rather 
large asymmetry between the two sides. Evidently, our 
results differ significantly from those of VM08 when 
considering the two sides of the stellar line separately. 
However, when merging signals from the two sides of the 
line, the corresponding absorptions seem comparable for 
the two analyses, yet this similarity (8.4% ± 1.8% here 
and in BJ07 compared to 7.3% ± 2.% in VM08 ) is only 
incidental. 

In a second step, we varied the wavelength domain 
of integration over the stellar line. As shown in Figure 
4, both our selected wavelength window and the one of 
VM03 give quite similar LCs with a transit absorption 
rate around ~ 8.4% ± 1.8% for our spectral domain and 
- 8.85% ± 2.2% (compared to 15% ± 4% in VM03 ) 
for the domain selected by VM03. In addition, when 
considering the whole stellar line, our absorption rate 
remains around ^ 8.9% ± 1.8%, consistent with our 
initial estimate and also with recent analysis using 
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Fig. 2. — (a) A sample spectrum of HD209458b Ly-a emission 
line observed by HST/STIS/G140M medium resolution after cor- 
rection for all sky background contaminations (thick line). The 
sky background contaminations (including the geocoronal and the 
interplanetary Ly-a emissions for the date of observation; see Ta- 
ble. 1) is displayed (thin line) to highligh its line wings contami- 
nation due to the effect of the line spread function of the instru- 
ment. The spectral domain of the sky background contamination 
is shown as a shaded area, (b) Sky and stellar integrated flux ver- 
sus distance (arcsec) displayed along the HST/STIS/G140M slit. 
The flux was evaluated by integration over the spectral domain 
([121.483, 121.643]) englobing the spectral extent of the sky emis- 
sion and part of the stellar line. The signal spike indicated by a 
shaded area around positions ~ 10" corresponds to the stellar line. 
The signal has been binned to eleven pixels corresponding to the 
default spatial extent of a point-like source. The signal gaps re- 
spectively around positions ~ — 1" and r^ 20" correspond to the 
location of the fiducial occulting bars of STIS. The variation (5%) 
of the signal along the slit that results from the non-uniformity 
of the STIS MAMA detector is comparable to the statistical noise 
level outside the spectral domain [121.541, 121.584] nm of contam- 
ination by the sky background signal. 

HST/ACS observations ()Ehrenreich et al.. I (|2008t l). 

What, then, is the physical origin of this weak sensitiv- 
ity of the absorption during transit to the selected wave- 
length range? Before we proceed to answer this question, 
in the following we study the time variability of the stel- 
lar in order to explore its potential impact on the transit 
absorption line profile. 

4. TIME VARIABILITY OF THE HD209458 Ly-o SIGNAL 
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Fig. 3. — Light curves obtained from our time series (as shown 
in Fig. 1 of BJ07) by integration of the signal in the spectral 
windows [121.483, 121.536] nm for the blue side (filled circle) of 
the stellar line and [121.589, 121.643] nm for the red side (open 
circle). Despite the noise level, the two curves behave similarly 
during transit with a total obscuration of ^ 8.5% ± 2.5% for the 
blue side and ~ 8.2% ± 2.5% for the red side (see text). If any 
difference, it would be in the opposite direction of VM08's claim, as 
indicated by the time bin around ~ 5000s where the red side seems 
more absorbed than the blue one. For reference, the beginning and 
end of the visible transit are shown by vertical dotted lines. 



-10 



Orbital phase angle (Degrees) 
-5 



1.20 



1.10 T 



3 1 .00 — 
q 

]c : 

g 

■S 0.90 - 



Y^ 



0.80 r ♦ BJL07 
: O VID03 

0.70 E \ . ^ 

-10000 



-5000 

Time from transit (s) 



5000 



Fig. 4. — Light curves obtained from our time series (as shown in 
Fig. 1 of B J07) by integration of the signal in the spectral windows 
[121.483, 121.643] nm (filled diamond, BJ07) compared to the one 
obtained from the wavelength windows [121.515, 121.61] nm (open 
diamond, VM03). A total obscuration of ~ 8.4% ± 1.8% is derived 
for the first domain and 8.8% ± 2.2% for the second one during 
the planetary transit. The closeness of the two absorption rates 
rules out any large spectral dependency as claimed by VM03 and 
VM08. For reference, the beginning and end of the visible transit 
are shown by vertical dotted lines. 

Time variability of the stellar signal is considered here 
as a potential cause for the discrepancy between the 
different data analysis related to the HST/STIS/G1401VI 
medium resolution observations. VM03 and VM08 
used the solar Ly-a variability as a proxy to catch the 
variability of the Ly-a emission line of HD209458 (solar 
type star). This approach may not be adequate because 
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Fig. 5. — HD209458 spectra obtained from data files 
(O6E203010, O6E203020, O6E203030) for difl'erent orbital po- 
sitions of the planetary companion. A "quiet" spectrum from 
the remaining HST data set (eg Table 1) is shown for refer- 
ence. The spectra are rebinned lay 2 at a spectral resolution of 
~ 26Kms~^. Several spectral features (indicated by arrows) ap- 
pear in the O6E2030x0 spectra but could not be found in the other 
exposures. Moreover, the strength of these features changed over 
time as if related to the orbital position of the exoplanet. However, 
because the phase positions are set by the HST observing schedule, 
the transit period is not fully covered, thus making it difHcult to 
correlate the signal variation with the transit effect and the related 
orbital positions of the exoplanet. 

no two stars of the same type behave similarly over short 
time scales. Moreover, HD209458 has a nearby orbiting 
exoplanet that may induce tidal and electromagnetic 
modifications of the stellar signal, effects that do not 
exist for our sun. As a matter of fact, HD209458 was 
previously suspected to have a relatively moderate chro- 
mospheric activity from Call H and K lines that were 
record ed over full orbits of the system (jShkolnik et al.. I 
(|2005( l). More recently, time variability of the stellar 
Ly-a flux was noted and showed both for a specific 
wavelength range and at different orbital phases (see Fig. 
2 & 3b in BJ07). However, the noise level was relatively 
high making difficult the diagnostic on the origin of the 
signal variability. Another possible way to characterize 
time variability is to carefully inspect the individual 
data set so far obtained by comparing the stellar signal 
at different phase positions to test any correlation with 
the orbital position of the planet. Such correlation 
may result from unkown electromagnetic interactions 
suspected to o perate between the pl anet and the star 
during transit (jCuntz et al.. I (|200Clf ): ICranmer fc Saar.1 
([2007)), much like what we know from the Jupite r-Io or 
Jupiter-Ganymede systems (jClarke et al.. I ()2002[ )). 

A careful inspection of the four visits thus far avail- 
able in the archive reveals that, indeed, the data sub- 
set obtained during the third transit of the VM03 ob- 
serving program (data files O6E203010, O6E203020, and 
O6E203030; see Table 1) show extra features that are 
not observed in the other visits of the HD209458 system 
(Fig. 5). Some of the observed features are well above 
the noise level and are suspected to correlate with the 
planet orbital position, particularly for the spectral fea- 
ture that appears around ^ 121.3 nm in the far blue side 
of the stellar Ly-a line. If we restrict the analysis to the 



solely Ly-a line profile in the range ±200 km/s from the 
line center, it is difficult to descriminate a time variabil- 
ity from statistical variations at low S/N . Yet we can 
estimate such time variability by comparing integrated 
fluxes from the different exposures as shown for example 
in Fig. 3 of VM03 and Fig. 2 of VM08. In these figures, 
the data points corresponding to the third transit (the 
so called C points) behave strangely with a departure of 
up to 3 (7 far from the theoretical transit light curve. 

All these indications tend to support that data set 
O6E2030X0 bears an unusual variability that is not 
found in the other exposures. The extra features are not 
related to any Earth's energetic particles activity along 
the HST orbit, as all observational parameters are com- 
parable to the other data sets. We suspect that the data 
set O6E2030x0 may have corrupted the diagnostic on 
the transit depth and absorption profile in VM03 anal- 
ysis. In our case, the impact of the same data set was 
largely diminished by the splitting of the exposures into 
our time series and by having extra observations (data 
files O4ZEA40X0, see Table 1) strengthening the weight 
(~ 25% more) of the pool of "quiet" exposures. Finally, 
the uniqueness of the disturbed data set O6E2030X0 and 
the noise level make it difficult to definitely assess a cor- 
relation with the orbital position of the planet, a ques- 
tion that remains open until future observation of the 
HD209458 system can be obtained after the FUV capa- 
bilities of the HST are restored. 

5. ABSORPTION LINE PROFILE DURING TRANSIT 

In practice, to obtain a self-consistent absorption 
profile, both steady-state unperturbed (out-of-transit) 
and in-transit line profiles are required and will be 
processed by iteration. In the present case, the iteration 
process converged in only a few steps, providing the 
absorption profile as shown in Figure 6. The obtained 
absorption spectral profile is typical of damped absorp- 
tion systems, though the spectral domain around the 
line center is hidden by, respectively, the interstellar 
medium absorption and the contamination from the 
sky background. The first thing to note is that the 
atmospheric absorption is statistically flat far into the 
stellar line wings down to 121.5 nm on the blue and 
up to 121.63 nm on the red, including the spectral 
position of the stellar line peaks. As the peaks' intensity 
provides the highest contribution to the LC shown in 
Figures 3-4, it becomes apparent that a flat absorption 
of - 8.9% ± 2.1% proposed in our initial study (BJ07) 
is a rather good approximation. Beyond that frequency 
domain, it is difficult to accurately assess the spectral 
variation due to the noise level. 

Binning the absorption line profile in wavelengths 
improves the signal quality, yet it does not reveal any 
asymmetry between the blue and red sides at a reduced 
spectral resolution of ^ 26km s~^. As previously shown 
in Figure 3, fiuxes integrated, respectively, on the blue 
and red sides of the emission line behave quite similarly 
over time, supporting the contention that both wings 
are absorbed almost symmetrically during transit in the 
spectral domain < 200km s^^ from line center. Any 
difference, if at all, must be weak and/or occur outside 
the available spectral domain such that current data 
noises and/or the sky background contamination hide 
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Fig. 6. — Absorption line profile (filled circles) during transit 
of HD209458b obtained as a ratio between Ly-a line profiles in 
transit and the unperturbed one. The sky background spectral 
domain as defined here is indicated by two vertical dashed lines, 
while the one of VM03 is indicated by a hashed area. The spectral 
window was restricted to ~ ±200kms^-'^ from line center, because 
the signal becomes much too noisy beyond. A fiat absorption rate 
of ~ 8.4% fits rather well the obta ined absorption profile. Using 
the standard atmospheric model of IMunoz. I 1 1200 7 1). a first model 
fit to the absorption line profile (solid line) has a Lorentzian shape 
with extended wings. Because data bear statistical and variability 
noises and the core spectral domain is absent at Earth's orbit, other 
scenarios of gas distribution may produce an absorption profile 
(dashed line) that fits rather well with the observed one, and thus 
cannot be ruled out. 
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In the following, we take a theoretical approach 
to verify whether such observational properties of 
the absorption profile are predictable, particularly 
its Lorentzian shape (Fig. 6). First, we recall some 
theoretical background on the opacity properties within 
a strong conservative line like Ly-a. Typically for a res- 
onance line, the absorption is dependent on wavelengths 
with strong opacities in the line core and much less in 
the (absorption) line wings. For a typical temperature 
of ~ lOOqOK f o r the HD209 458b inflated atmosphere 
(jMunoz. I ()2007D : lYelle7l(|2004f )). the absorption Doppler 
width is almost ^ 0.0052 nm (or ^ 13km s^^), which is 
comparable to the spectral resolution (~ 13km s^^) of 
the HST observations used here. Referring back to the 
HD209458 line profile (VM03, BJ07), the fine's peaks 
appear at ^ ±100kms~^, or equivalently, ^ ±8 Doppler 
units far from the line center, placing them well in the 
wings of the absorption line profile of the planetary 
atmosphere. Because of the steep fall off of the flux in 
the stellar wings, the contribution of the signal around 
the peaks' positions is dominant in any estimate of the 
flux attenuation during transit. Also, one should note 
that for a target that is not spatially resolved, such as 
transiting exoplanets, the observed signal extinction is 
a weighted average of absorption from different regions 
of the exoplanetary extended atmosphere, with the 
outer layers having the maximal weight. This could be 
examined by writing the basic definition of the observed 
(transmitted) stellar line profile during full transit: 



Fiv) = exp{-T^p,e,v))F,„iv,e).p.dp.d9 (!) 

where $(p, u)) is the absorption profile by H atoms 
versus the velocity v (or wavelength) from the line 
center, Fin{v) is the unperturbed stellar line profile, and 
t{p, 9) is the opacity along the line of sight that goes 
through the cloud's volume at the impact parameter p 
from the star's disc center and a polar angle 6. The 
integration is evaluated over the whole extent of the 
stellar disc, allowing various distributions of the gas 
around the planet that may lead to the same observed 
extinction. This degeneracy would not appear if the 
atmospheric layer were thin compared to the size of the 
planet. 

To illustrate the multiplicity of solu tions, a s tandard 
atmospheric model of HD209458b (jMunoz. I (J2007); 
lYelle. I (|2004l )) was used. As shown in Figure 6, a rather 
good fit to the absorption line p rofile dur i ng tra nsit is 
obtained for the DIVl model of lMunoz.l ()2007[ ). The 
DIVl model includes heavy elements contribution to 
the H production and appears up-to-date regarding 
chemis t ry rea ctions. If using the standard model of 
lYelle. I ()2004D that does not include heavy elements 
chemistry, the same good fit is obtained when the H 
abundance is scale d by a constant facto r Beyo r id the 
difference between iMunoz. I ()2007l ) and lYelle. I (|2004D 
models, the scaling constant may also reflect other 
processes of H atomic production that are not accounted 
for in classical photochemistry models. In the past, 
the study of the thermosphere of giant planets such 
as Jupiter and Saturn and their satellites showed the 
need for such enhancement of the H production in 
order to fit the Ly-a observations thus far obtained by 
Voyager 1 fc 2, lUE, Ga fi leo, o r HST (see, for i n stance , 
Ben-Jaffel et"aL~l (l2007l [T995h : iClarke et al.TI (fT99lh : 
Sandel et al.. I (|1980f )). For example, in the case of 
Jupiter, a planetwide enhancement (^ 3 times) of the H 
production is expected from magnetospheric protons and 
electrons precipitation in the polar regions (mainly 112 
or heavy molecules dissociation by energetic particles). 
For reference, the same process produces the classical 
auroral emissions. The constituents produced by auroral 
processes and the input energy and momentum in the 
polar regions are then tra nsported and distribute d 
over all latitudes /long itudes (jSommeria et al.. 1 (|1995| ): 
iBougher et al.. I ()2005D ). As of today, the H abundance 
and distribution o n Jupi ter are not yet fully uncovered 
(jBen-Jaffel et al.. I ()2007| )). The idea here is that the 
magnetospheric energy input in the auroral regions may 
play a key role in the upper atmosphere of HD209458b's 
composition and structure. Such effects are not taken 
into account in atmospheric models thus far published 
in the literature. On that basis, for a planet such as 
IID209458b orbiting at 0.05AU from its parent star, 
thereby receiving much more energy through its polar 
regions than expected on Jupiter, the conditions may 
favor a rather strong enhancement of the H production 
that may lead to larger atmosph e ric op acities at Ly-a 
than shown in Fig. 15 of lMunoz 
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In the model calculations shown in Figure 6, a Voigt 
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function was assumed for the H atomic absorption profile 
in order to account for the natural broadening of the Ly- 
a line. As expected from the line shape, the atmosphere 
is opaque and consiste nt with the line -of-sight opacities 
used (see Figure 15 in MunozTI ()2007l )). The line shape 
reflects the effect of the broadening in Ly-a line wings 
by the large H column along the line of sight that leads 
to a large and relatively flat absorption far from the line 
center. However, it is important to point out that the in- 
terpretation of the absorption level versus wavelengths is 
not straightforward because the absorption is a weighted 
contribution from different atmospheric layers located at 
different distance from the planet's center (e.g.. Equation 
1). For that reason, we stress that the fit to the absorp- 
tion profile and the corresponding atmospheric opacities 
correspond to one possible sketch of the H distribution 
that may not contain all of the key physical processes at 
play in the real spatially extended system. Other sce- 
narios of gas distribution may lead to the same observed 
absorption profile (e.g. Fig. 6). For example, we are 
able to find a more extended atmosphere for which the 
extinction by the exoplanet atmosphere fits rather well 
with the observed absorption profile (see Fig. 6). As 
it is impossible to obtain any information on the line 
core because of the interstellar medium absorption and 
the sky background contamination, only very high S/N 
FUV observation at a level better than (0.1 — 1.0)% may 
disentangle the different processes in play. 

6. COMPARISON WITH PREVIOUS ANALYSIS 

At this point, we are able to compare our result with 
previous analysis that used most of the same data set 
(VM03, VM08). The immediate question is, what could 
be the origin of the difference in the transit absorption 
rate as derived for the same wavelength domain DWl : 
[vl,v2] = [121.515, 121.61] nm-namely, 8.8% ± 2.2% 
in this study-compared to 15.% ± 4% in the case of 
VM03 and VM08? VM08 claimed that a difference 
in the selected wavelength domains was the cause of 
the difference in the absorption rates. To support 
their claim, those authors derived a transit absorption 
of 15% ± 4% from their analysis for the wavelength 
domain DWl, while a different transit absorption of 
7.3% ± 2.1% was derived for the wavelength domain 
DW2 : [v3,v4] = [121.484, 121.643] nm. The problem is 
that the coherency claimed by VM08 is merely illusory 
because the provided rates are simply inconsistent. 
Indeed, one may note that the considered wavelength 
window DWl is embedded in DW2, with DWi : 
[vl,v3] & [v4,v2] = [121.484, 121.515] nm & [121.61, 
121.643] nm as the difference domain. Considering 
the out-of-transit stellar line profile as provided by 
VM03 and VM08, one may wonder: how could an 
absorption rate possibly be cut by a factor of two when 
extending the wavelength domain from DWl to DW2 
and only adding the falling-off signal from the very 
limited DW3 stellar wings? In other words, what could 
be the transit absorption rate in the complementary do- 
main DW3 if 15% is obtained in the first domain DWl 
and 7.2% in the second one DW2 as estimated by VM08? 

Starting from the definition of the transmitted fiux in 
Equation 1, we write that the absorption Ajjwi inte- 
grated over the wavelength domain DWl is a weighted 



sum of the absorption A]j]Y2 integrated in the domain 
DW2 and Adw3 in the domain DW3. For a sake of 
clarity we introduce the following notation for the signal 
attenuation A^w corresponding to a prescribed wave- 
length domain DW : 



A 
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WW 
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is the fraction of transmitted signal with 
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In above equations, the subscripts N stands for numera- 
tor and D for denominator, and Fin is the unperturbed 
stellar fiux. Note that the contribution to the integra- 
tion from the sky background wavelength domain is set 
to zero. 

Starting with the most extended DW2 domain, we may 
write: 



I 



DW2 



J^^F,„{v)dv 



(5) 



Splitting the integration in the numerator over the dif- 
ferent domains, we obtain 



J„3 F,n[v)dv 



(6) 



Taking into account the definition of the rates in the 
DW2 and DWi domains, we can write: 



LDW2 



Idwi X Idwi + Idw3 X I] 



DW3 



^DW2 



(7) 



We now may derive the absorption in the DW3 wave- 
length domain : 



Idw3 



1dw2 X Ijj-W2 ^ Idwi x I^wi 



^DW3 



(8) 



All integrals ijji^ are known and directly derived from 
the unperturbed (out-of-transit) line profile. Now, us- 
ing the numbers provided by VM03 and VM08 for the 
attenuation in the different wavelength domains of inter- 
est here, namely that Adwi = 15% ± 4% (or Idwi = 
85% ± 4%) and Adw2 = 7.3% (or Idwi = 92.7% ± 2%), 
we can derive the fraction of the signal transmitted in 
the DW3 wavelength domain should be Idw3 ~ 108%. 
This means that in order to make consistent the 15% and 
7.2% rates derived by VM08 respectively for the DWl 
and DW2 spectral domains, an enhancement of 8% of 
the stellar signal is required during transit in the wave- 
length window DWS. Surprisingly, when the absorp- 
tion rate proposed by VM03 is divided by a factor 1.08 
(8% enhancement), then we recover an absorption rate of 
8.7% for the DWl spectral domain, consistent with our 
present estimation (see Table 2). We believe that the 
normalization used by VM03 to compensate stellar vari- 
ability by using the stellar wings signal as a proxy may be 
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TABLE 2 

Sensitivity of the transit absorption rate to the spectral domain 



Wavelength window (nm); reference 
[121.515,121.610]"; {VM03) 
[121.483, 121.643]''; (BJ07) 
[121.483,121.536]; (BJ07, blue side) 
[121.589,121.643]; (BJ07, red side) 
[121.515,121.550]; (VM03, blue side) 
[121.580, 121.610]; (VM03, red side) 



Published absorption (rate, reference) 

15% ± 4% {VM03) 

8.9% ± 2.1% (BJ07); 7.3% ± 2.1% (VMOS) 

9.8% ± 1.8% (VM08) 

5.2% ± 1.0% (VM08) 



Absorption rate (estimated here) 
8.8% ± 2.2% 
8.4% ± 1.8% 
8.5% ± 2.5% 
8.2% ± 2.5% 
9.3% ± 3.3% 
8.3% ± 2.9% 



one of the causes of the discrepancy between their results 
and ours. Indeed, VM03 and VMOS repeatedly claimed 
that because most of the absorption should, in their es- 
timation, occur in the line core, then the stellar Ly-a 
line wings (out-domain in VMOS definitions) should not 
show any spectral absorption. As shown here (e.g., Fig. 
6), the transit absorption profile has very extended wings 
with a substantial absorption up to the out-domain used 
by VMOS, which casts doubt upon their approach to es- 
timating the transit attenuation and to handling stellar 
time variability using the line wings' signal as a reference. 

7. DISCUSSION 

Despite serious limitations of current technology to 
resolve exoplanets from their parent stars, we learned a 
lot on their atmospheric structure thanks to the transit 
technique. For the extrasolar planet HD209458b, we 
are close to sketching its atmospheric structure over 
several scales with detection, respecti vely, of neutral 
sodium in the lower atmosphere (C harbonneau et al.. I 
(|2002| )). very extended corona of hydrogen in the upper 
atmos phere, and an in t ermed iate layer of temperature 
jump (jBallester et al.. I ()2006f )). For exoplanets orbiting 
very close to their host stars, these results confirm the 
scenario of atmospheric inflation by the huge stellar 
flux, a fact that raises the question of atmospheric 
evaporation and the survival of these planets over 
several stages of stellar evolution. 

Therefore, it may be worthwhile to discuss the 
adequacy of thus-far proposed models to explain the 
hydrogen distribution around hot giant exoplanets in 
the light of the constraints derived here and to propose 
new hints that may help investigate this fundamental 
problem. Since VMOS disseminated their results on 
HD209458, many studies followed regarding the mass 
loss from exoplanets, particularly that class of planets 
that closely orbit to their parent stars. Initially, VMOS 
proposed evaporation followed by radiation pressure 
acceleration of hydrogen atoms beyond the Roche lobe 
leading to a comet-like shape of the hydrogen nebula 
that escape away from the star. As we know from BJ07 
and the present study, this scenario is not confirmed by 
the HST/STIS Ly-a observations because no signature 
of asymmetry could be found as initially claimed by 
VMOS. The se cond difficulty of the VM OS model was 
pointed out bv iHolmstrom et al.. I (|2008f ). namely, that 
the radiation pressure has not enough time to operate 
before the atoms are ionized. 

These conclusions do not mean that evaporation 
is not operating in the HD209458 system. In fact, 
the only process being dismissed is the acceleration 



by radiation pressure because, consistently with the 
conclusions of the present analysis, it does not seem to 
be efficient enough to produce the necessary population 
that must leave a cle ar spectral asymmetry i n the Ly-a 
absorption profile. ISchneiter et al.. I ()2007D tried SD 
hydrodynamical simulation of the hydrogen nebula 
around IID209458 and predicted several light curves 
for the planetary transit, depending on the escape 
efficiency of the atmosphere. Unfortunately, these 
authors underestimated the absorption by the hydrogen 
cloud by assuming a Gaussian absorption profile for the 
H atoms at Ly-a , an assumption that is not adequate 
for this atomic transition. In all cases, the natural 
width of the upper level of the H atom induces a natural 
broadening, described by a Lorentz profile t hat must 
be included in th e abso r p tion calculations (Mi halasTI 
(0,978); Ben-Jafl[el et al., I ©88)). When included, this 
natural broadening should significantly enhance the 
absorption in the Ly-a line wings and consequently the 
absorption for the total line. If the velocity distribution 
of atoms is Maxwelian, then the attenuation is enhanced 

— 2 

by the ratio ~ —^ — r of a Voigt to a Gaussian function 
that may be several order of mag nitude beyond x ~ S 
Doppler units from the line center (jNeufeld. I (119901 )). 



In a later study, IHolmstrom et al.."1 (|2008D proposed 
that the inflated atmosphere of IID209458b should be 
extended beyond the magnetopause region of the planet 
so that charge exchange with the stellar wind protons 
would necessarily produce energetic neutral atoms 
(ENAs). In this process, the size of magnet o spher ic 
cavity (obstacle region in IHolmstrom et alTI ()2008l )) 
is one of the key parameters that defines the ENAs 
content. The more the obstacle region is extended, 
the less ENAs are produced. Posing the magnetopause 
position at ~ 4.2Rp, these authors were able to account 
for the HD209458b transit's absorption features as 
claimed by VMOS assuming a stellar wind with a density 
of 2 X lO'^cm"'^ and a bulk velocity ^ 50 kms~^. As 
far as the conclusions of VMO S and VMOS hav e been 
invahded, the fit proposed in IHolmstrom et al.. I ()2008D 
is no more adequate and one may question if other 
parameters of the stellar wind could account for the 
constraints derived in BJ07 and in the present study. 

Before we conclude on this question, it is important to 
discuss some dif ficulties of the ENAs prod uction model 
as proposed by IHolmstrom et al.."1 ()2008D . Because of 
the nature of the charge exchange reaction, we must 
have two populations: one is the ENAs and the second 
is a family of ions flowing outward in the same direction 
as the planetary atoms that produced them. In their 
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Fig. 7. — (Dashed line) Ve l ocity field of neutral hydrogen pro- 
posed by "Holms trom et al.. I II2008I ) as a best fit to account for 
ENAs production through charge exchange between HD209458b 
atmospheric hydrogen and impinging stellar wind protons. (Solid 
line) The same distribution when accounting for natural broaden- 
ing due to the finite life time of the H second atomic level. The 
net eff'ect of natural broadening is to spread the initial distribution 
over a wider velocity range and making its red side stronger (~ 5 
times) and almost symmetric compared to its blue side in the ref- 
erence frame of the star. This result is independent of the intensity 
of the distribution that should be normalized to unity to derive the 
local atomic absorption profile. 



publication, iHolmstrom et al.."1 (|2008f ) just removed 
that population from their simulation box. Also, the 
stellar wind should slow down when approaching the 
assumed obstacle. A self-consistent formulation of the 
problem should allow for a cascade of ENAs and ions 
productions through t he ch arge exchange process (see 
lOsterbart and Fahr. I (|1992f) ). Such ions will induce 
their own ram pressure against the stellar wind, pushing 
away the magnetospheric cavity (the obstacle) out- 
ward at larger distances (Ra tkiewicz and Ben-Jaffel. I 
(|2002l )). How great that distance is depends on self- 
co nsistent simulations that n eed to be done. According 
to IHolmstrom et al.. I (|2008l ) (Fig. 7 in their supple- 
mentary draft), moving the obstacle (magnetopause) 
further out will reduce the ENAs production and, 
consequently, the attenuation in the red and blue 
part of the spectrum. The conclusion here is that the 
current version of the the ENAs model as proposed by 
IHolmstrom et aLA (2008) lack the coherent treatment 
of magnetospheric structure and response that is much 
needed to evaluate whether this process is important 
in shaping the hydrogen distribution around HD209458b. 



Another difhculty the IHolmstrom et al.. I (|2008f ) model 
faces is the formulation of the attenuation process of 
the Ly-g emissi o n due the atomic hydrogen. As in 
ISchneiter et al.. I ()2007f ). these authors evaluated their 
attenuation profile during transit while overlooking the 
fact that the second level of a hydrogen atom has a 
finite life time that induces a natural broadening of the 
Ly-a line. The natural broadening is represented by a 
Lorentz function that describes the absorption profile of 



H atoms in their rest frame. Therefore, for a gas with 
a gi ven velocity distribu t ion (l ike the one shown in Fig. 
2 of 'Holm strom et al.. I (|2008f )). the absorption profile 
that should be used for attenuation calculations must be 
obtained as the convolution be tween the Lorentz profile 
and the velocity distribution (jMihalas. I (1978) ). For 
the s ake of clarity, we r eprod uced the velocity profile 
from IHolmstrom et al.. I (|2008( ) and convolved it with a 
single Lorentz profile for a gas temperature of 8000K. In 
general, when dealing with a non-homogeneous medium, 
for every position across the planetary disc (see Equation 
1), one needs to evaluate the atomic absorption profile 
locally for the corresponding temperature. The same 
caution is required when dealing with several popu- 
lations of particles. As sho wn in Fig. 7, the velocity 
distribution field derived bv IHolmstrom et al.. I (|2008) ■ 
when compared with the same distribution including 
the effect of natural broadening, leads to distinct results 
with stronger and extended line wings and a smoothing 
out of the red side, making the new distribution much 
symmetric in the star's reference frame than initially set. 
Finally, we remark that the effect of natural broadening 
is even more visible when the ENAs contribution is 
smaller. 

With that in mind, one may well imagine a hybrid so- 
lution where a thin layer of superthermal atoms (ENAs 
production, turbulence energized particles, etc.) is on 
top of a classical thermal atmosphere that may absorb 
the Ly-a stellar flux in the way it is observed (e.g.. Fig. 
6). The [H] thickness of the thin layer should not ex- 
ceed few percent of the column of the underlying atmo- 
sphere. This configuration takes its origin in similar sce- 
narios that also occur in the so -called bulge region in th e 
upper atmosphere of Jupiter ("Ben-Jaffel et al.n (|1993f) : 
Hommcria et al., (1995): Emcrich et al.. I (|1996l )). As 
previously note d, however, ENAs and i ons production 
as proposed bv IHolmstrom et al.. I (|2008f ) must be eval- 
uated in conjunction with the evolving size of the mag- 
netopause cavity so that the effective ENAs population 
that contributes to the Ly-a absorption is consistent with 
the magnetopause position and the gas structure around. 
Also, the corresponding absorption profile should ac- 
count for the effect of natural broadening of the Ly-a line 
as described here (eg Fig. 7). In practice, many solutions 
may be found that may require spanning the stellar and 
atmospheric models parameters in a self-consistent way. 
This effort is out of the scope of this study. 

8. CONCLUSION 

HST archive observations of the Ly-a emission of 
HD209458 are revisited to test the absorption rate sen- 
sitivity to the spectral range considered inside the Ly-a 
emission line and to derive the absorption line profile 
during transit. Our first conclusion is that the absorp- 
tion depth during transit is weakly dependent on the 
assumed spectral range assumed, with rates in the range 
of (8.4 — 8.9)% ± 2.% for all possible wavelength windows 
as reported in our initial study (BJ07). Self-consistent 
derivation of the absorption line profile during transit 
shows a symmetrical Lorentzian shape, typical of op- 
tically thick conditions of the hydrogen nebula around 
HD209458b. Whether using direct light curves from the 
time series or inspecting the absorption line profile dur- 
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ing transit, neither method shows any indication of pre- 
ferred absorption on the blue side of the stellar line in 
the spectral domain r^ ±200kms~^ from line center. Fi- 
nally, we stress that because the planet is not resolved 
spatially, several scenarios of the gas distribution may 
lead to the observed absorption line profile. In that con- 
text, two scenarios that fit rather well to the absorption 
profile are suggested. A first model assumes an optically 
thick atmosph ere of HD20945 8b obtained from the stan- 
dard model of lMunoz, ! (2007) . A second class of models 
consists of a standard model that may have a thin layer 
of superthermal neutrals (ENAs, turbulence energized 
particles, etc.) on top of it, much like what we observe 
for th e upper atmosphere of Jupiter (jBen-Jaffel et al..l 
(|1993f) : lEmerich et al.. I (|1996D V In all cases, the natu- 
ral broadening that was missed in most past studies is 
required in the modelling of the Ly-a absorption profile. 
Moreover, for Enas production, a self-consistent treat- 



ment of the charge exchange process in conjunction with 
the size of the magnetospheric cavity is much needed to 
test the validity of such process. In any case, high resolu- 
tion observations with a S/Nhetter than (0.1 — 1.0)% are 
required to efficiently constrain the hydrogen distribution 
and its time and phase variations along the planetary or- 
bit. 
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